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Abstract 

We report results on the decay B~ — >• DcpK~ and its charge conjugate using a data sample 
of 85.4 million BB pairs recorded at the T(4S') resonance with the Belle detector at the KEKB 
asymmetric e^e~ storage ring. Ratios of branching fractions of Cabibbo-suppressed to Cabibbo- 
favorcd processes are determined to be B{B- -> D^K-)/B{B- D^tt') = 0.077 ± 0.005(stat) ± 
Omeisys), B{B- DiK~)/B{B~ Dnr-) = 0.093 ± 0.018{stat) ± 0.008(sys) and B{B- 
D2K-)/B{B- D2-K-) = 0.108 ± 0.019(stot) ± 0.007(sys) where the indices 1 and 2 represent 
the CP=+1 and CP=— 1 eigenstates of the — system, respectively. We find the partial- 
rate charge asymmetries for B~ — > DcpK~ to be Ai = 0.06 it 0.19{stat) it 0.04(sys) and A2 = 
-0.19 ± 0.17{stat) ± 0.05(sys). 

PACS numbers: 13.25.Hw, 14.40.Nd 
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The extraction of 03 [ij, an angle of the Kobayashi-Maskawa triangle 0, is a challenging 
measurement even with modern high luminosity B factories. Recent theoretical work on B 
meson dynamics has demonstrated the direct accessibility of 03 using the process B~ 
DK~ 0, 0|- If the is reconstructed as a CP eigenstate, the b ^ c and b ^ u processes 
interfere. This interference leads to direct CP violation as well as a characteristic pattern 
of branching fractions. However, the branching fractions for D meson decay modes to CP 
eigenstates are only of order 1 %. Since CP violation through interference is expected to be 
small, a large number of B decays is needed to extract 03. Assuming the absence of — D^ 
mixing, the observables sensitive to CP violation that are used to extract the angle 03 
are, 
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where the ratios R^^''^ and R^° are defined as 
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Di and D2 are CP-even and CP-odd eigenstates of the neutral D meson, r = \A{B- — >■ 
D^R-) /A{B- D^K-) \ is the ratio of the amplitudes of the two tree diagrams shown in 
Fig. 1 and 6 is their strong-phase difference. The ratio r corresponds to the magnitude of 
CP asymmetry and is suppressed to the level of ~ 0.1 due to the CKM factor (~ 0.4) and 
a color suppression factor (~ 0.25). Note that the asymmetries Ai and A2 have opposite 
signs. 




K- 



B- 




K- 



FIG. 1: B- D^K- and D^K~. 

The ratio of the Cabibbo-suppressed decay B^ —>■ D^K~ to the Cabibbo-favored decay 
B- D'^n- has been reported by CLEO to be R^' = 0.099l[]:E]lt t^oe while Belle finds 
rd'' = 0.079 ± 0.009 ± 0.006 0. Assuming factorization, the ratio R^° is expected to be 
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ta.v? 6c{fKl I-kY ~ 0.074 in the tree-level approximation, where 6c is the Cabibbo angle, 
and fx and /tt are meson decay constants. The measurements are in good agreement with 
this theoretical expectation. 

Previously, Belle reported the observation of the decays DiK^ and D2K~ 

with 29.1 fb~^ 0. This paper reports more precise measurements of these decays with a data 
sample of 78 fb~^, containing 85.4 million BB pairs, collected with the Belle detector at the 
KEKB asymmetric-energy e~^e~ (3.5 on 8 GeV) collider operating at the T(4S') resonance. 
At KEKB, the T(45') is produced with a Lorentz boost of = 0.425 nearly along the 
electron beamline. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a three- 
layer Silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
Silica aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised 
of CsI(Tl) crystals located inside a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the coil is instrumented to detect 
mesons and to identify muons (KLM). The detector is described in detail elsewhere |9j]. 

We reconstruct D° mesons in the following decay channels. For the flavor specific mode 
(denoted by Df), we use D° K't:^ [lOj. For CP =-|-l modes, we use Di K^K^ and 
7r~7r+ while for CP =-1 modes, we use D2 K^g'^^, K^g(j), K^guj, K^rf and K^i . 

The charged track, and 7r° selection requirements have been described in Ref . Q • For 
each charged track, information from the ACC, TOF and specific ionization measurements 
from the CDC are used to determine a K/tt likelihood ratio P{K/tt) = /{Lk + L^^), where 
Lx and L,^ are kaon and pion likelihoods. For kaons (pions) from the K^tt^ mode we 

used the particle identification requirement of P{K/'k) > 0.4 (< 0.7). For kaons from the 
— ^ K~K^ mode we require P{K/tt) > 0.7 while for pions from tt^tt^ mode we 

require P{K/n) < 0.7. 

The u mesons are reconstructed from tt+tt^tt'^ combinations in the mass window 
0.732 GeV/c^ < M(7r+7r"7r°) < 0.82 GeV/c^ with the charged pion particle identification 
requirement P{K/7r) < 0.8. To reduce the contribution from the non-resonant background, 
a helicity angle cut |cos ^heil > 0.4 is applied where 6'hei is the angle between the normal to 
the uj decay plane in the u rest frame and the u momentum in the rest frame. To remove 
the contribution from —>■ K*~p^, we require the Kg7T~ invariant mass to be greater than 
75 MeV/c^ from the K*~ nominal mass. 

The mesons are reconstructed from two oppositely charged kaons in the mass window 
of 1.008 GeV/c^ < M{K+K-) < 1.032 GeV/c^ with P{K/n) > 0.2. We also apply the 
hehcity angle cut |cos 6*1101! > 0.4 where 6'hei is the angle between one of the daughters 
in the rest frame and the momentum in the rest frame. We form candidate rj 
and rj' mesons using the 77 and rjn^n' decay modes with mass ranges of 0.495 GeV/c^ < 
M(77) < 0.578 GeV/c2 and 0.903 GeV/c^ < M(r/7r+7r-) < 1.002 GeV/c^, respectively. 
The 7] momentum is required to be greater than 0.5 GeV/c. Both rj and 77' candidates 
are kinematically constrained to their nominal masses. The candidates are required to 
have masses within ±2.5(T of their nominal masses, where a is the measured mass resolution 
which ranges from 4.9 MeV to 17.7 MeV depending on the decay channel. A mass and 
(wherever possible) vertex constrained fit is then performed on the remaining candidates. 

We combine the and it~/K~ candidates (denoted by h) to form B candidates. We 
apply tighter particle identification cuts, P{K/tt) > 0.8(< 0.8) for prompt kaons (pions), 
to identify B~ D^K~{tt~) events. The signal is identified by two kinematic variables 
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calculated in the center-of-mass (cm.) frame. The first is the beam-energy constrained 

mass, Mbc = v-Ebeam ~ \Pd where po and ph are the momenta of and K~ /ti~ 

candidates and -Ebeam is the beam energy in the cm. frame. The second is the energy 
difference, /S.E = Ed + — -Ebeam, where E^, is the energy of the candidate. Eh is 
the energy of the /n^ candidate calculated from the measured momentum and assuming 

the pion mass. Eh = \J\ph\'^ + '^n^- With this definition, real D^n^ events peak at 

AE = even when they are misidentified as D^K~ , while B^ D^K^ events peak 

around AE = —49 MeV. Event candidates are accepted if they have 5.2 GeV/c^ < Mbc < 
5.3 GeV/c^ and \AE\ < 0.2 GeV. In case of multiple candidates from a single event, we 
choose the best candidate on the basis of a determined from the differences between the 
measured and nominal values of Md and Mbc. 

To suppress the large combinatorial background from the two-jet like e~^e~ qq {q = u, 
d, s or c) continuum processes, variables that characterize the event topology are used. We 
construct a Fisher discriminant F, from 6 modified Fox Wolfram moments Furthermore, 
cos 9b, the angle of the B flight direction with respect to the beam axis is also used to 
distinguish signal from continuum background. We combine these two independent variables, 
F and cos 6b to make a single likelihood ratio variable (LR) that distinguishes signal from 
continuum background. We apply a different requirement for each sub-mode based on the 
expected signal yield and the backgrounds in the Mbc sideband data. For B~ D^n' 
where — * K~tt^, K^K^ we require LR > 0.4 whereas for tt+tt^, K^tc^, Kg(j), 

KgUJ, Kg-T] and K^t]' we require LR > 0.6. To give an example of the performance of this 
selection, the LR > 0.4 requirement keeps 87.5 % of the B~ D^[—>- i^'~7r"'"]7r~signal while 
removing 73 % of the continuum background. 

The signal yields are extracted from a fit to the AE distribution in the region 
5.27 GeV/c^ < Mbc < 5.29 GeV/c^. The B^ D^ti~ signal is parameterized as a double 
Gaussian with peak position and width fioated. On the other hand, we calibrate the shape 
parameters of the B~ D^K~ signal using the B~ D^7^~ data. This accounts for the 
kinematical shifts and smearing of the AE peaks caused by the incorrect mass assignments 
of prompt hadrons. The peak position and width of the B^ D^K~ signal events are 
determined by fitting the B" D^tt~ distribution using the kaon mass hypothesis for the 
prompt pion, where the relative peak position is reversed with respect to the origin. The 
shape parameters for the feed-across from B^ D^tt~ are fixed by the fit results of the 
B^ D^7i~ enriched sample. The continuum background is modeled as a first order poly- 
nomial function with parameters determined from the AE distribution for the events in 
the sideband region 5.2GeV/c^ < Mbc < 5.26GeV/c^. Backgrounds from other B decays 
including contributions from B~ — > D*^K~ and B^ —>■ D^K*^ are modeled as a smoothed 
histogram from Monte Carlo simulation. The fit results are shown in Fig. 2. 

The ratios of branchin g fr actions of Cabibbo-suppressed to Cabibbo-favored processes 
are determined as follows |l2j |. 

^ N{B- DK-) ^ r^jB- ^ Dn-) ^ ^ 
N{B- Dn-) r]{B- DR-) e{K) ' 

where N is the number of observed events, rj and e are the signal detection and the prompt 
pion/kaon identification efficiencies, respectively. The signal detection efficiencies were deter- 
mined from a Monte Carlo simulation e.g. ri{B~ — > DfTT~) = 44.6 % and ri{B~ — >■ DfK~) = 
42.5 %. The particle identification efficiencies for the prompt pion and kaon, e(7r) and ^{K), 
are determined from a kinematically selected sample of D*^ D^[—>- /^^7r+]7r+ decays, 
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TABLE I: Signal yields, feed-acrosses and ratios of branching fractions. The errors on R are 
statistical and systematic, respectively. 



Mode 




B Dtt- 
events 


B- DK- 
events 


B Dtt- 
feed-across 


B{B-^D0K-) 
^ ~ B(B-^D'^n-) 


B- 
B- 
B- 


Dfh- 
Dih- 
D2h- 


6052 ± 88 
683.4 ± 32.8 
648.3 ± 31.0 


347.5 ± 21 

47.3 ± 8.9 

52.4 ± 9.0 


134.4 ± 14.7 
15.6 ± 6.4 
6.3 ± 5.0 


0.077 ± 0.005 ± 0.006 
0.093 ± 0.018 ± 0.008 
0.108 ± 0.019 ± 0.007 


TABLE IL Yields, partial-rate charge asymmetries and 90 % C.L intervals for asymmetries. 


Mode 




N{B+) 


N{B-) 


Acv 


90 % C.L 


B^^ 
B^^ 


DfK^ 
DiK^ 

D2K^ 


165.4 ± 14.5 
22.1 ± 6.1 
29.9 ± 6.5 


179.6 ± 15 
25.0 ± 6.5 
20.5 ± 5.6 


0.04 ± 0.06±0.03 
0.06 ± 0.19 ±0.04 
-0.19 ± 0.17±0.05 


-0.07< Af <0.15 
-0.26< Ai <0.38 
-0.47< A2 <0.11 



where the K~ and tt"'" mesons from D° candidates have been selected in the same cm. mo- 
mentum (2.1 GeV/c < Pcm. < 2.5 GeV/c) and polar angle regions as prompt hadrons in the 
B~ — > Dh~ decay. With our requirement P{K/t:) > 0.8, the efficiencies were determined 
to be e{K) = 0.768 ± 0.001 and e(7r) = 0.976 ± 0.001, and the rate for misidentification of vr 
as K is 0.024 ± 0.001. The ratio of B ^ Dn' feed-across to B ^ Dn- signal is 2-2.5 % 
which is consistent with the measured pion fake rate. The ratios of Cabibbo-suppressed to 
Cabibbo-favored decay modes are shown in Table L The double ratios are found to be 

7^l = 1.21 ± 0.25{stat) ± 0.14(sys), 
7^2 = 1.41 ± Q.27{stat) ± 0.15{sys) 

for CP-even and CP-odd eigenstates, respectively. The systematic errors in the ratios 
are due to the uncertainty in yield extraction (3-7 %) and particle identification (1 %). 
The systematic error in the yield extraction includes uncertainties in the BB background 
and signal shape parametrization. The uncertainty in the AE signal shape parametrization 
was determined by varying the mean and width of the double Gaussian parameters within 
their errors. The uncertainty from the slope of the background was determined by changing 
its value by its error. Both of the resulting changes were included in the systematic error 
from fitting. Also other backgrounds including rare decays such as B~ — > K~K^K^ and 
B" — > ir~7r+7r~, which could contribute to the AE signal region, are estimated from the 
sideband data. This uncertainty (0.5-3 %) is also included as a source of systematic error. 

The asymmetries Ai^2 are evaluated using signal yields obtained from separate fits to the 
B^ and B^ samples shown in Fig. 3. The results are given in Table II. We find 

Ai = 0.06±0.19(stai) ±0.04(sys), 
A2 = -0.19 ± Q.n{stat) ± 0.05(sys) 

where the systematic uncertainty is from the intrinsic detector charge asymmetry (3.2 %), 
the B~ and S+ yield extractions (2.4-3.3 %), and the asymmetry in particle identification 
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FIG. 2: AS distributions for (a) DfK-, {h)B- DfK~, {c)B~ Z^ivr", (d)-B~ DiK' , 

{e)B~ — > D2TT~ and {i)B^ — > D2K^ . Points with error bars are the data and the sohd lines show 
the fit results. 

efficiency of prompt kaons (1 %). The intrinsic detector charge asymmetry is calculated from 
the B~ —>■ Z)°[— »• K~7t~^]tt~ sample. The systematic error from yield extraction is calculated 
by changing the fitting parameters by ±la. 

In summary, using 78 fb~^ of data collected with the Belle detector, we report measure- 
ments of the decays DcpK~, where Dqp are the neutral D meson CP eigenstates. 
These supersede the results reported in 0,13. The ratios of the branching fractions R^^'^ for 
the decays B~ DcpK~ and B~ — >■ DcpT^~ are consistent with those for the fiavor specific 
decay within errors. The measured partial-rate charge asymmetries Ai^ are consistent with 
zero. 
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FIG. 3: AE distributions for the charge conjugate modes (a) B DiK , {h)B^ DiK^ , 
{c)B- D2K-, {d)B+ D2K+. 
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